A graphics processing unit-based robust numerical model for solute transport driven by torrential flow condition

Solute transport simulations are important in water pollution events. This paper introduces a finite volume Godunov-type model for solving 4×4 matrix form of the hyperbolic conservation laws consisting 2D shallow equations and equations. The adopts Harten-Lax-van Leer-contact (HLLC)-approximate Riemann solution to calculate cell interface fluxes. It can deal well with changes dry wet interfaces an actual complex terrain, it has strong shock-wave capturing ability. Using monotonic upstream-centred scheme (MUSCL) linear reconstruction slope Runge-Kutta time integration method achieve second-order accuracy. At same time, introduction graphics processing unit (GPU)-accelerated computing technology greatly increases speed. is validated against multiple benchmarks, results good agreement analytical solutions other published numerical predictions. third test case uses GPU central (CPU) calculation models which take 3.865 s 13.865 s, respectively, indicating that increase speed by 3.6 times. In fourth case, comparing calculated traditional CPU, efficiencies under different resolution grids 9.8–44.6 times higher than those CPU. Therefore, better potential previous large-scale simulation solute incidents. provide reliable theoretical basis data support rapid assessment early warning accidents.